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Abstract 



We have recently presented evidence that in configurations dominating the regular- 
ized pure-glue QCD path integral, the topological charge density constructed from the 
' overlap Dirac operator organizes into an ordered space-time structure. It was pointed 

out that, among other properties, this structure exhibits two important features: it 
is low-dimensional and geometrically global, i.e. consisting of connected sign-coherent 
p) ' regions with local dimensions 1 < d < 4, and spreading over arbitrarily large space- 

^ • time distances. Here we show that the space-time structure responsible for the origin 

of topological susceptibility indeed exhibits global behavior. In particular, we show 
numerically that topological fiuctuations are not saturated by localized concentrations 
^ ■ of most intense topological charge density. To the contrary, the susceptibility saturates 

only after the space-time regions with most intense fields are included, such that geo- 
metrically global structure is already formed. We demonstrate this result both at the 
fundamental level (full topological density) and at low energy (effective density). The 
drastic mismatch between the point of fluctuation saturation (~ 50% of space-time 
at low energy) and that of global structure formation (< 4% of space-time at low en- 
ergy) indicates that the ordered space-time structure in topological charge is inherently 
global and that topological charge fluctuations in QCD cannot be understood in terms 
of individual localized pieces. Description in terms of global brane-like objects should 
be sought instead. 

1. Introduction. Understanding the nature of topological charge fiuctuations is an essential 
step toward clarifying the structure of the QCD vacuum. Indeed, not only are there impor- 
tant physical phenomena associated with topological charge (such as rj' mass, 6'-dependence 
and possibly spontaneous chiral symmetry breaking (SChSB)), but there are also intriguing 
theoretical reasons to believe that topological charge structure in QCD may provide a partic- 
ularly clear window on the gauge theory/string theory connection, since topological charge 
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is dual to Ramond-Ramond charge in type IIA string theory [1] . Nevertheless, it is only re- 
cently that a first-principles direct study of space-time topological structure in a well-defined 
lattice-regularized setting was carried out [2, 3]. This development was made possible by 
important advances in understanding lattice chiral symmetry [4, 5, 6] which, among other 
things, led to the construction of novel lattice topological charge density operators [7, 5]. 
These operators are not just convenient in that the corresponding topological susceptibility 
can be computed without any subjective manipulations of the gauge field and that theoretical 
properties of the lattice topological field match those in the continuum [8]. In addition, they 
exhibit features that are inherent to smooth backgrounds on continuum manifolds, but quite 
unexpected for general backgrounds contributing to the regularized QCD path integral. In 
particular, the index theorem is exactly satisfied relative to the lattice Dirac operator used 
to define the density [7], and the global charge of a generic configuration is strictly stable 
with respect to small local changes of the gauge field. 

In previous work we have studied the space-time behavior of such topological density 
based on the overlap Dirac operator in pure-glue QCD [2]. We found a clear excess of or- 
dered structure compared to random distributions of charge. In fact, topological charge 
organizes into two oppositely-charged low-dimensional "sheets" that fill a macroscopic frac- 
tion of space-time, and arc built around a connected "skeleton" of points with high intensity ^. 
This structure is embedded in space-time in a non-trivial geometric manner and it should be 
regarded as fundamental since it includes fiuctuations at all scales up to the lattice cutoff. ^ 
We note that a viable candidate for the fundamental structure in typical configurations con- 
tributing to the QCD path integral must reproduce the correct behavior of the topological 
density correlator [3]. A particularly notable requirement is that it has to conform to the 
fact that {q{0)q{x)) < at arbitrary non-zero distance [10]. Indeed, for the above two- 
sheet/skeleton structure the non-positivity of the correlator at non-zero distances arises in 
an ordered manner due to the fact that layers of low-dimensional structure with alternating 
sign permeate the space-time. 

Two basic aspects of the topological charge structure described in [2] arc particularly 
intriguing and possibly far-reaching. First is the strictly low- dimensional nature of the struc- 
ture ^. By that we mean that the structure only contains parts with local dimensions strictly 
less than four ^. This suggests some interesting possibilities for the way in which gauge fields 
guide the propagation of pions, and will hopefully shed new insight into the relation between 
topological charge fluctuations and SChSB. The second aspect is the geometrically global 
nature of the structure. This expresses a geometric fact that the subset of space-time oc- 
cupied by the structure (sign-coherent sheet or skeleton) is a connected set spreading over 
maximal available distances. In this paper we focus on the distribution of charge within the 
geometric structure and show that the global behavior is necessarily present in the physically 

^From the geometric point of view, both the sheets and the skeleton behave somewhat similarly to Peano's 
space- filling curve [2]. 

^Analogous structure was recently reported to be seen also in the 2-d CP(N-l) models [9]. 

^The low-dimensional nature of fundamental topological field is reflected to some degree also in low-lying 
Dirac modes [11, 12]. More precise form of this correspondence is yet to be understood. The notion of 
strictly low- dimensional structure is also emerging using indirect projection techniques [13]. Its relation (if 
any) to the structure in topological field is yet to be understood. 

"^The dimension is at least one since the sign-coherent sheet/skeleton is a connected global structure. 
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relevant substructure responsible for the origin of topological susceptibility. In particular, 
we demonstrate that the bulk of topological susceptibility in QCD cannot be carried by 
localized concentrations of most intense topological charge. Just to the contrary, we will 
argue quantitatively that the space-time structure carrying topological charge fluctuations is 
inherently global (super-long-distance), in the sense discussed in Ref. [14]. In other words, our 
results indicate that this structure (both at the fundamental level and at low energy) cannot 
be broken into localized objects that could be used as a basis for valid physical analysis of 
the origin of topological susceptibility in QCD. Analysis in terms of global low-dimensional 
branc-likc objects should be sought instead. 

2. The Argument. The methodology for analysis of global behavior in space-time structure 
of local observables has been developed in Ref. [14]. Let us first review the basic ingredients of 
this approach in case of topological charge density. In what follows, we will assume that q{x) 
is a well-defined lattice topological density in the sense that the topological susceptibility can 
be directly computed via the associated correlator or, equivalently, via x = HQ^) ~ 
where Q = J2x If D is the overlap Dirac operator [6] based on the Wilson-Dirac kernel 
with mass —p, then the density [7] 

q{x) = ^tr-f^D^^^ = -tr75(l - ^^^,3,) (1) 

satisfies this requirement. 

The basic idea of Ref. [14] is to simultaneously monitor the physics (in this case topological 
susceptibility) and the geometry (in this case global behavior) as the space-time structure in 
q(x) is built starting from the most intense point and then gradually lowering the intensity 
threshold for points included in the structure. Formally, we define the sets S'^{f) (fraction 
supports) containing a fraction / = N{S'^{f ))/N{Q) of the most intense points as ranked by 
\q{x) I . Here N{T) denotes the number of points contained in the subset F of discretized torus 
Q. The fraction / thus plays the role of the monitoring variable. We emphasize that what 
wc mean by "structure" here is the partition of iS^(/) into maximal path-connected ^ subsets 
Sf{f) containing the points with topological density of the same sign [14, 2]. This definition of 
geometric structure emphasizes connectedness because extended connected concentrations of 
strong fields can naturally facilitate the long-distance physics (e.g. lead to inherently global 
Dirac eigenmodes [14] and Goldstone boson propagation). Considering sign coherence is 
physically motivated by the fact that quarks of preferred chirality are attracted to regions 
of definite sign. 

To monitor the saturation of topological fiuctuations, we compute the cumulative function 
C{f) of topological susceptibihty [14], which represents the fraction of the total susceptibility 
carried by a fraction / of the most intense points, i.e. 

c(/) . |W ,(/) . mm-iQinr 

where 

Qif) = Yl ?(^) V = «'^(^) (3) 

xesi{f) 

^By path-connected lattice set we mean a set of points any pair of which can be connected by a nearest- 
neighbor path passing only through points within the set. 
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are the charge associated with the fraction support and the total volume respectively. To 
monitor the global behavior in the structure, wc compute two associated characteristic func- 
tions. The first is the maximal linear size imaxif) of the connected component, i.e. 

Cax(/) = max{£(5f(/))} (4) 

Here the linear size i{r) for arbitrary F C f2 is the maximal distance between two points of 
the set, i.e. i{T) = max{ \x — y\ : x,y E T}. The second characteristic is the point-wise 
average of the linear size over the points of the structure. More precisely, to x e Sf{f) we 
assign the linear size if{x) = i{Sf{f)) of the corresponding connected component, and take 
the average over 

The motivation for introducing imaxif) ^^ptaif) is that their values signal different regimes 
of global behavior in the structure. In particular, when imaxif) reaches values imaxif) ~ ^(^) 
(close to the "size of the space-time") at some typical fraction fmax, then fmax characterizes 
the onset of global behavior. In other words, an intensity threshold that corresponds to 
/ ^ fmax is low cuough to iucludc not only the most intense local concentrations of the 
charge distribution, but also the connecting ridges which join them into a global structure. 
In addition, for fractions / > fmax: where the largest connected piece is global, we would 
also like to know at what value of / most of the points in the structure are part of the global 
component. The function iptaif) characterizes this property and its satTiration at some value 
fpta > fmax indicates that global behavior is prevalent throughout the structure. 

We computed the above characteristics for the ensemble (51 independent configurations) 
of Iwasaki gauge action [15] on an 8^ lattice with scale a = 0.165 fm determined from string 
tension. The physical volume of the system is thus V — 3 fm^. The topological charge 
density operator (1) with p = 1.368 (k = 0.19) has been used to compute the site-by-site 
values of full topological charge density. The details of numerical implementation for overlap 
matrix-vector operation can be found in Ref. [16]. The results are shown in Fig. 1. 

The data indicates that imaxif) ^ strict (full) saturation point at fmax ~ 0.14. Also, 
'iptaif) exhibits a clear saturation "knee" at fpta ~ 0.20. At the same time, the cumulative 
function of susceptibility saturates more gradually and reaches a full saturation point (within 
the statistics) at fc ~ 0.80. The central message of Fig. 1 is that at fractions where the 
geometric structure consists of localized pieces (i.e. iptaif) — imaxif) <^ 1) ^ only a small 
percentage of the total topological susceptibility is accounted for. In other words, topological 
susceptibility is not dominated by locahzed carriers of topological charge. Moreover, the 
susceptibility saturates only at fractions obviously larger than those at which the structure 
becomes prevalently global (i.e. / > fpta)- This result quantitatively justifies the proposition 
put forward in Ref. [2], that the physically relevant vacuum structure in topological charge 
density is inherently global (super-long- distance). In other words, our data indicate that it is 
not correct to view vacuum fluctuations of topological charge in terms of localized lumps of 
topological charge, while it strongly supports the idea that global extended objects should 
be used basis for valid physical analysis. 

It should be noted here that our arguments started with the assumption that space-time 
regions containing the most intense fields (in this case topological charge density) are most 

^In what follows we will always quote distances in units of £{^l) , the maximal possible distance. 
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V = 8^ Iwasaki action a = 0.165 fm 51 configs 
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Figure 1: The functions C(/) (monitoring topological susceptibility) and imaxif), ^ptaif) 
(monitoring the degree of global behavior) for full topological density in ensemble with renor- 
malization group improved Iwasaki gauge action. Distances are measured in units of 

relevant for the physics in question (in this case topological susceptibility). This assumption 
is a posteriori justified by the fact that the computed C(/) is a concave function on its 
domain within the errorbars. In other words, the ordering of space-time points by intensity 
results in monotonically decreasing contribution to susceptibility as the fraction is increased. 



3. Low Energy. Our analysis in the previous section was carried out for full topological 
charge density, i.e. with fluctuations at all scales up to the lattice cutoff included. An at- 
tractive feature of topological charge density operators based on Ginsparg- Wilson fermionic 
actions is the possibility of defining an effective topological field [3] via the eigenmode expan- 
sion of Eq. (1), i.e. 

A 



(6) 



A<A^a 



where c'^(x) = ^(x)75'?/''^(x) is the local chirality of the mode with eigenvalue A. The 
effective field (x) describes topological fluctuations up to energy scale defined by the 
probing chiral fermion. Since zero modes are included in the expansion for arbitrary A^, we 
have J2xQ{^) = Q = J2xQ^ = — n+. Here ^^-(+) is the number of zero modes with 
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negative (positive) chirality. This has three noteworthy imphcations [3]. (a) The topological 
susceptibility is independent of A^. (b) Similarly to q{x), the effective density is a perfect 
topological field in the sense that (for generic backgrounds) its global sum is strictly stable 
with respect to small local changes of the gauge field, (c) (x) satisfies the index theorem 
on the lattice in the same sense that q{x) does [7]. 

The introduction of the effective topological density opens the possibility of describing the 
space-time structure of topological fluctuations in a scale-dependent manner as is appropriate 
in field theory [3]. The emerging precise relation between the structure at low energy and 
the fundamental structure will be described elsewhere. Here we focus on the issue of global 
behavior at low energy and how it compares to the behavior at the fundamental level. To 
do so, we have constructed effective densities at = 600 MeV for three ensembles of 
Wilson gauge action with parameters described in Table 1. Lattices at different cutoffs have 
the same physical volumes V ^ 3 fm^. The quoted values of lattice spacing were obtained 
via the Sommer parameter using the interpolation formula given in Ref. [17]. For each 
configuration in a given ensemble we have calculated the zero modes and a fixed number of 
complex-conjugate eigenmode pairs specified in Table 1. The Ritz variational method [18] 
was used for this calculation. A typical accuracy of calculated eigenvalues (as measured by 
differences of complex-conjugate pairs) was one part in 10^ or better. The largest cutoff at 
which the effective density could be constructed for £3 was just above A^ = 600 MeV (larger 
for £1 and £^2)- 

We have calculated the cumulative function of topological susceptibihty C{f) and the 
global characteristics imaxif) and iptaif) for the three ensembles with the results shown in 
Fig. 2 (left column). One can see that the main change at low energy compared to the 
fundamental structure is that the functions imaxif) iptaif) saturate even faster. For 
example, imaxif) ^3 = 0.055 fm) reaches full saturation point at fmax < 0.06. We thus 
conclude that the inherently global nature of the topological structure in the QCD vacuum 
is even more apparent at low energy. Moreover, the wide range of cutoffs spanned by our 
data shows convincingly that this conclusion will not change in the continuum limit. 

To see the relation between the degree of saturation in susceptibility and in global be- 
havior more directly, we follow Ref. [14] and eliminate the "monitoring variable" / in favor 
of C. This defines the functions imaxi^) and iptaiC) that quantify the degree of global be- 
havior for any given value of susceptibility saturation. These functions arc plotted in the 
right column of Fig. 2 and exhibit the strongly concave behavior characteristic of inherently 
global structure. From the data one can see that the structure becomes prevalently global 
when it carries only C ~ 0.2 of total susceptibility. 

To quantify the vast discrepancy between the fractions of space-time volume characteriz- 
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Table 1: Ensembles of Wilson gauge configurations. Zero modes and fixed number of lowest 
complex eigen-pairs were computed for each configuration. 
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Figure 2: The cumulative function of topological susceptibility and the global characteristics 
of the geometric structure at low energy. See the discussion in the text. Distances are 
measured in units of 
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Wilson action V = 3 fm^ = 600 MeV 
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Figure 3: The continuum extrapolation of /°'^ associated with imaxif), ^ptaif) and 
The hnearly extrapolated values are quoted in the plot. 



ing the appearance of global behavior and the saturation of susceptibility, we have calculated 
the scale dependence of fractions at which £max{f), £pta{f) and C(/) reach the reference sat- 
uration value y{f^'^) = 0.9. In Fig. 3 we plot f^'^{a) for the three characteristic functions. In 
all three cases we performed a linear extrapolation to the continuum limit with the extrapo- 
lated results quoted in the plot. This calculation shows that 90% saturation of susceptibility 
requires inclusion of about 50% of space-time occupied by most intense topological charge 
density in the continuum limit. At the same time, the 90% saturation of £max{f) and iptaif) 
requires less than 2% and 4% of space-time respectively. This represents a rather striking 
demonstration that topological structure in QCD vacuum is inherently global (super-long- 
distance) at low energy. 

4. Conclusions. The goal of the study described here was to justify the proposition 
of Ref. [2] that the ordered vacuum structure in topological charge density, first observed 
there, is based on underlying g'/o 6a/ connected object(s) present in the vacuum. The original 
suggestion of Ref. [2] was based on the purely geometrical observation that the regions of 
space-time containing the most intense topological density have a strong preference to or- 
ganize into a global long-range structure via the formation of a skeleton '', i.e. the minimal 

^In the language developed in Ref. [14] the formation of a skeleton corresponds to saturation of ipta- This 
definition is simpler than the one used in Ref. [2] and leads to similar quantitative results. 
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hard-core substructure exhibiting the prevalently global behavior. In order to ascertain that 
a description of physics entails the existence of global geometric objects rather than some 
hypothetical localized objects forming the skeleton, it is necessary to evaluate the contri- 
bution of such objects to the relevant physical observable. Ref. [14] developed a formalism 
for performing such a calculation in the case of topological susceptibility. We adopted this 
formalism and the results presented here demonstrate that the global behavior dominates 
at decisively smaller fractions of space-time than those that are necessary to saturate the 
topological susceptibility. This is true both at the fundamental level and at low energy where 
the difference of relevant fractions is strikingly large. From these observations we conclude 
that the vacuum topological structure behaves as inherently global (super-long- distance) in 
QCD. ^ The relevance of this result resides in the fact that it directly indicates that a vahd 
theoretical description of topological charge fluctuations in QCD vacuum needs to be based 
on global brane-like objects. 
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